hasil analisis di bab menunjukkan bahwa C, (T,j) dan CT, j) yang diberikan oleh rumus 
(7.7) dan (7.8) menjadi tak-terdefinisikan bila titik T berhimpit dengan PG) atau P-1). 
Dengan pengertian tersebut kita membuat anggapan bahwa panel ke j cukup diwakili oleh 
titik tengahnya, bila dipandang dari segi hitungan distribusi kecepatan di sepanjang 
permukaan aerofoil. Bila T berada di tengah-tengah panel PO(j) maka khusus untuk panel 
PO) rumus (7.7) dan (7.8) disederhanakan menjadi (lihat lampiran bab 5). 

C,(T, j) —0 (9.15) 

C,(T,j) 0 (9.16) 
Sekarang marilah kita bahas metoda panel yang berdasarkan pada pemodelan dengan panel- 
panel vorteks dengan distribusi kekuatan seragam dan kondisi batas Dirichlet. Untuk kasus 
ini nilai fungsi arus di titik T diberikan oleh rumus berikut 


wWT)- YT. cosa- XT.sinat- 31 )C,(TI) (7.17) 


ja 
dimana C, (T,j) adalah koefisien pengaruh panel PO) pada nilai fungsi arus di titik T, yaitu 


CT, j) - PO(i) tORIT, j).In OTT, j)— PR(T,).In PT(T,j) - TR(T, j). ACT, j) 
(7.18) 
Kondisi batas Dirichlet yang untuk metoda panel sumber diberikan oleh persamaan (7.11), 
untuk kasus metoda panel vorteks diubah menjadi 
JP--1 


»Cki,i).x() — Di) untuk 1 JP (7.19) 


j4 
dimana 
(1,4) - C,(ij) 


xi) ti) 


| untuk 1xixs JP 
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D(i) - 2x(XP/i).sina — YP(i).cosa) 
Cli, JP 41) —-—2x 
KIP 41) Wae 
Apabila semua nilai y(j) sudah diketahui maka distribusi kecepatan di sepanjang permukaan 


aerofoil dapat dihitung seperti telah dijelaskan untuk metoda panel sumber, yaitu 


V(T) - Iv (1) “op (9.20) 


x 


dimana 


1 ANTA 
VA(T) - cosa t— Y'C,(T, j).y(i) 
21 45 


1 Apa 
V,(T) - cosa | ON CAT, j).y(i) 
Tj 


C,(T,j) --A(T,j).cos B(j) ca SA Jain B(5) 





Gk) RI Jeon) Atm) 


Perhatikan bahwa prosedur di atas itu cukup rumit dan melibatkan banyak perhitungan. 
Khusus untuk metoda panel vorteks, kita dapat mengurangi upaya komputasi yang 
dibutuhkan untuk menghitung V(T) dengan menerapkan pengertian tentang kekuatan 
vorteks dan sirkulasi. 

Bayangkan sebuah elemen segi panjang yang digambar mengurung suatu elemen 
panjang As pada panel PO) (lihat gambar (7.1)). Sisi AB dari segi panjang ABCD berada 
sedikit di bawah dan sejajar dengan panel PO(j), sedangkan sisi CD berada sedikit di atas 
dan sejajar panel PO(j). Sisi AD dan BC tegak lurus pada PO(j) dengan panjang mendekati 
nol. Kecepatan aliran sepanjang sisi DC adalah V sedangkan di sepanjang sisi AB kecepatan 
bernilai nol karena AB berada di dalam bentuk aerofoil. Kecepatan di sepanjang sisi AD dan 
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BC bernilai nol karena kedua sisi tersebut normal pada PO(j). Sirkulasi di sepanjang kurva 
ABCD adalah 
Tancp “ Vas: AB 4 Voc. BC— V,... CD — V,,,.DA 

Tanep “04#0— V.CD—0 ——V.CD ——V.As 
Kekuatan vorteks pada elemen panjang As adalah yAs dan karena kekuatan vorteks itu 
sama dengan sirkulasi maka 

Txn»cp #—V.As— y.As 
atau 

V--—y (7.21) 
Pembahasan di atas berlaku untuk elemen panjang pada permukaan atas aerofoil. Untuk 
permukaan bawah aerofoil hasil yang diperoleh adalah 

V-—y (7.22) 
Seperti dapat dilihat, untuk kasus metoda panel vorteks kita tidak perlu menghitung nilai 
kecepatan dengan rumus (7.20) karena persamaan (7.21) dan (7.22) menunjukkan bahwa 
kecepatan di sepanjang panel PO(j), yaitu Vj), ternyata sama dengan kekuatan vorteks per 
satuan panjang, y(j). Dengan pengertian teersebut sekarang kita dapat menerapkan kondisi 
Kutta untuk metoda panel vorteks sebagai berikut. Kondisi Kutta menuntut bahwa 
kecepatan di tepi buritan permukaan atas itu harus sama dengan kecepatan pada permukaan 
bawah. Kecepatan pada panel PO(1) di permukaan bawah adalah V(1)- sedangkan pada 
panel PO(JP) di permukaan atas kecepatan adalah V — —y(JP). Karena V(d) harus sama 
dengan V(JP) maka kondisi Kutta adalah 

1G) --y(JP)-—0 (7.23) 
Persamaan (7.23) dapat ditulis menjadi bentuk berikut (lihat persamaan (7.12)) 

JP 


HC IP 41, j).x(j) — DJP «-1) (9.24) 


ja 
dimana 


C(JP #1, j) —0 untuk 24 j£ JP-1 danj- JP-1 
C(JP 1,1) — C(JP 4-1, JP) —1 
D(JP 4-1) -0 
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Persamaan (7.19) kemudian digabung dengan (7.24) menjadi 


JP 


5 cli,i).x(j) Di) untuk 1Ki« JP (7.25) 


Seperti dapat dilihat metoda panel vorteks itu jauh lebih sederhana dan lebih efisien 
dibandingkan dengan metoda panel sumber. Disamping itu metoda panel vorteks dapat 
digunakan untuk memodelkan aliran di sekeliling aerofoil bergaya angkat. Sekarang marilah 
kita pelajari masalah pemodelan metoda panel vorteks dengan kondisi batas Neumann, 
dimana nilai turunan fungsi dan bukan nilai fungsi itu sendiri yang ditentukan pada batasan 
medan solusi. 

Untuk masalah kondisi batas Neumann kita akan menerapkan persyaratan bahwa udara 
tidak boleh menembus keluar masuk permukaan aerofoil dengan menuntut bahwa komponen 
kecepatan arah normal pada permukaan aerofoil bernilai nol. Mengingat bahwa kecepatan 
normal adalah 


-.& 
Wa (7.26) 


maka nilai dari turunan fungsi potensial ditentukan pada batasan medan solusi atau 
permukaan aerofoil dan oleh karena itu permukaan masalah disebut masalah syarat batas 
Neumann. Arah normal pada panel P@(i) adalah 


li) — —sinBli)x4-cosP(i)y (9.27) 
jadi nilai komponen kecepatan normal pada panel P@(i), di titik T(i) adalah 
V, (T(i)) - V,(T(i)) cosB(i) — v, (T(i))simp(i) (7.28) 


Persamaan di atas dapat diuraikan lebih lanjut menjadi 


V 


n 


(i) —sina. cosB(i) — cosa.sin Bli) 
"elestisi)-cosgti) —e,ti,i)-simp(G(i) 
Apabila T(i) dipilih berada di tengah-tengah panel PO(i) maka 


(7.29) 


erli,j) 
PT(i,j) 





CM4)-AG) co) Jamata untuk i-j 
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C,li,i)-0 


C,li,j) — ag Ta 


C,li,i) 0 





Jon )-A(i,4).sinB()) untuk izj 


“ 
& 
kami 
Ts 
— 


Perhatikan bahwa rumus untuk C,(i,i) dan C,(i,i) di atas berlaku untuk kasus dimana T(i) 


benar-benar berada pada panel PO(i). Apabila titik T(i) sedikit di atas PO(i), yaitu pada jarak 
O" dari permukaan aerofoil PO(i) maka 


C,li,i) --x.cosgli) 
C,li,i) —-x.singli) 


Kondisi batas Neumann sekarang dapat ditulis menjadi 


Cli,j).x()) —D(i) untuk 1XiS JP (9.30) 


Ms 


1 


dimana untuk 1 Xi « JP 
Clij) — C,(i, j).cosBli) —C,(i, i).sin Bli) 
x(4) —y(3) 


Dj) — 2x cos O&.sin B(i) — Sina.cos B(i)) 


Disamping harus memenuhi kondisi batas Neumann, model masalah juga harus memenuhi 
kondisi Kutta berikut 


11) #y(JP) —0 (7.31) 
yang dapat ditulis ulang menjadi bentuk 


Cli,)).xG) —Dli) untuk i- JP 1 (9.32) 


Ms 


J 


dimana 
C(JP -41,j) -0 untuk 25 js JP-1 
C(JP 1,1) - C(JP 4-1, JP) —1 
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D(JP 4-1)-0 
Persamaan (7.30) dan (7.32) dapat digabung menjadi 


Cl j).x() —D(i) untuk 1 Kis JP-H1 (7.33) 


ind 


U 


3 
Perhatikan bahwa sistem persamaan (7.33) terdiri dari JPt1 persamaan dan hanya 
melibatkan JP variabel anu, yaitu x(j). Ini berarti bahwa sistem persamaan (7.33) adalah 
sistem dengan ketentuan berlebihan (overdetermined system of eguations). Untuk membuat 
sistem persamaan tersebut menjadi sistem yang imbang, maka kita harus mengurangi jumlah 
persamaan dengan membuang salah satu dari persamaan (7.30), yaitu mengubahnya menjadi 
berlaku hanya untuk 1SIX JP-1. Selanjutnya kondisi Kutta (7.31) kemudian yang 
setimbang adalah 


Saida) -Dli) untuk 1Kix JP (9.34) 
ji 


dimana untuk 1 Si X JP-1 dan 14 j« JP 
Cli,)) — C,li,j).cosBli)—c,(i,j).sin Bli) 
xi) -yli) 
D(j) — 2nicosa. sin Bli) —sina. cosg(i)) 
sedangkan untuk i-JP 
C(JP,j) —0 untuk 2 js JP-1 
C(JP,1) - C(JP, JP) —1 
D(JP)-0 
Suatu pilihan lain yang dapat dilakukan untuk mengatasi masalah sistem persamaan 
berlebihan adalah dengan menambah jumlah variabel anu. Disamping panel-panel permukaan 


aerofoil, kita tentukan juga bahwa model masalah juga melibatkan satu panel vorteks 
tambahan di bisektor sudut ekor dengan panjang sembarang tetapi cukup kecil, misalnya 


0.002 atau 0.02. Panel ini disebut panel PO(JP#1) dengan kekuatan y(JP-4-1), yang 


merupakan variabel anu tambahan. Kondisi Kutta tidak dinyatakan daam bentuk persamaan 
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(7.31) tetapi dengan menuntut bahwa kecepatan aliran di titik tengah panel POUP-1), 
normal pada panel, harus bernilai nol. Koordinat titik O(JP41) atau P(JP#2) dapat dihitung 
dengan menggunakan hasil pembahasan sebelumnya (lihat persamaan (7.12)), yaitu 





0 ..al AXOP(1).AYOP(JP) -AYOP(1).AXOP/(JP) 
PER Se ana takar arora) 


(ap 41) — (pi) #p(3P)) 


XP(JP 4-2) — XP(JP 4-1)--0.01 


YP (JP --2) — tanB(JP 4-1).0.014 YP(JP 4-1) 


Sistem persamaan gabungan syarat batas Neumann dan kondisi Kutta, yang harus 
diselesaikan adalah 


3 di,i).xi) -Dli) untuk ISIS JP-H1 (9.35) 
j-l 


dimana untuk 1Si£ JP dan 1 j« JP-#1 

Ci, i) - C,i, i).cosB(i) —C,(i,i).sin li) 

xi) —yli) 

D(j) — 2nlcosa. sin B(i) —sin a.cosP(i)) 
Dengan pendekatan yang sama dengan yang telah dijelaskan di atas, metoda panel 
berdasarkan model dengan panel doblet dapat juga diturunkan, tetapi tidak akan dibahas di 
sini. Sekarang marilah kita bahas masalah metoda panel berderajat lebih tinggi, yaitu orde 
2, yang melibatkan distribusi kekuatan linear sepotong-sepotong. Dari pembahasan 
sebelumnya kita tahu bahwa kekuatan vorteks per satuan panjang itu ternyata sama dengan 
kecepatan aliran di sepanjang panel di titik yang sama (lihat persamaan (7.21) dan (7.22)). 
Dalam metoda panel derajat 1, kekuatan per satuan panjang dari singgularitas yang 
digunakan dianggap terdistribusi secara seragam di masing-masing panel. Untuk metoda 
panel vorteks, anggapan distribusi tersebut berarti bahwa kecepatan aliran bernilai tetap di 
sepanjang masing-masing panel dan terjadi loncatan mendadak dalam nilai kecepatan aliran 
dari satu panel ke panel lainnya. Dalam kenyataannya distribusi kecepatan aliran di sepanjang 
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permukaan aerofoil itu bersifat mulus, tanpa patahan-patahan di setiap titik ujung panel (lihat 
gambar (7.2)). Distribusi kecepatan yang mulus tersebut jelas dapat ditaksir dengan lebih 
baik dan tepat oleh distribusi linear sepotong-sepotong di masing-masing panel. Inilah logika 
yang mendasari pengembangan metoda panel derajat lebih tinggi. Konsep metoda panel 
vorteks derajat 2 dengan kondisi batas Dirichlet dapat dijelaskan sebagai berikut. 

Nilai fungsi arus di titik T yang diimbas oleh panel vorteks PO dengan distribusi kekuatan 
linear adalah 


W(T,20) --—(CWT,P0).YP 4 CAT,P0).(1O-YP)) (239 


dimana 


CU T,PO) - PO --OR.InOT—PR.InPT—TR.A 


"0 2 2 2 2 
POH2.PR PO'—PR'4#TR OT PO'— TR Imp7-—PR-TR 
4 2.PO 2.PO PO 
Untuk titik T(i) yang diwakili oleh indeks i dan panel PO(j) yang diwakili oleh indeks j, 
persamaan (7.36) dapat ditulis sebagai 


C2(T,PO) A 








visi) — ti). 126) C26,i)-(roti1) ye) (137) 
yang dapat ditulis ulang menjadi 

w(i,i) seleri,). tcCOfi,i).yo(i 1) (7.38) 
dimana 


CP(i,j) - cili,j)-cz(i,j) 


COfi,j) - czli,i) 
Kondisi batas Dirichlet untuk kasus distribusi kekuatan linear sekarang dapat ditulis sebagai 
berikut 


LN) co) Pn 


DT Woro 5 2x( XP(i).sin @a— YP(i).cosa) 
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Persamaan di atas dapat ditulis ringkas menjadi 


Sat, ).x) —Dli) untuk 1KiX JP (3.35) 
Na | 
dimana 
Cli,1) - CPli,1) 
Cli,j) - CP(i,j)) -COfi,j-1) untuk 2 j« JP 
4-1) - coli, JP) 


x(j) -yP(j) untuk 1S j JP-41 
X JP #2) Wae 
Di) - —x/(XP(i).sin A— XP(i).cosa) 
Kondisi Kutta untuk kasus ini dapat ditulis sebagai berikut 
YP(1) #YP(JP 4-1) —0 
Tetapi perlu diingat bahwa titik P(1) dan titik P(JP#1) itu berhimpit untuk aerofoil berekor 
lancip, jadi 
YP(1) — YP(JP 41) 
Kedua persamaan di atas hanya dapat dipenuhi apabila 
YP(1) — YP(JP 4-1) —0 (7.41) 
Pemasukan kondisi Kutta (7.41) ke persamaan (7.40) dapat memberikan hasil berikut untuk 
Isis JP 


JP 
5 cli,).x()) -Cli, JP 4-2).x(JP 1-2) — Dli) (9.42) 
i-2 
Persamaan (7.42) dapat ditulis lebih ringkas menjadi 
JP 
S Cli,)).x() —D(i) untuk 1 Ki « JP (7.43) 
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apabila ditentukan bahwa 

Cli,j) - CP(i,j) -CO(i,j-1) untuk 2 j« JP 

Cli,1) ——27 

x()) —YP(j) untuk 2« js JP 

XL) - Waero 

D(i) — —21(XP(i).sina — YP(i).coso) 
Perlu diperhatikan bahwa X(1) dalam persamaan (7.43) bukan sama dengan yP(1) dan nilai 
YP(1) tidak perlu dihitung karena kondisi Kutta (7.41) menentukan bahwa yP(1)-0. Sistem 
persamaan (7.43) adalah sistem yang setimbang karena terdiri dari JP persamaan dan 
melibatkan JP variabel anu, yaitu nilai w,.,, dan YP(j) untuk 2 j « JP Seandainya aerofoil 
berekor tumpul maka titik batas Dirichlet (7.40) dapat diterapkan untuk 14 js JP41. 
Selanjutnya, kondisi Kutta tidak diberikan oleh (7.41) tetapi oleh 

YP(1) -YP(JP 1) —0 
yang dapat ditulis ulang menjadi 

Jpa2 

S clii).x(i) —D(i) untuki — JP42 

ia 
dimana 

Cli,j))-0 untuk 2 jX JP danj- JP-42 

Cli,1) — Cli, JP--1) —1 

D(i)-0 
Ini berarti bahwa untuk aerofoil tumpul, sistem persamaan setimbang yang harus diselesaikan 
adalah 


Sati, H)-x() —D(i) untuk TKI JP 42 (9.45) 


ja 
dimana untuk 1S i« JPt2 


Cli,1) - CP(i,1) 
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apabila ditentukan bahwa 

Cli,j) — CP(i,j) -CO(i,j-1) untuk 2 j« JP 

Cli,1) -—27 

x)— yP(j) untuk 2S j« JP 

X(1) - Waero 

Di) - —2n(XP(i).sin a— YP/i).cosa) 
Perlu diperhatikan bahwa X(1) dalam persamaan (7.43) bukan sama dengan yP(1) dan nilai 
YP(1) tidak perlu dihitung karena kondisi Kutta (7.41) menentukan bahwa yP(1)-—0. Sistem 
persamaan (7.43) adalah sistem yang setimbang karena terdiri dari JP persamaan dan 
melibatkan JP variabel anu, yaitu nilai w,,,, dan YP(j) untuk 2 £ j « JP Seandainya aerofoil 
berekor tumpul maka titik batas Dirichlet (7.40) dapat diterapkan untuk 14X j £ JP 1. 
Selanjutnya, kondisi Kutta tidak diberikan oleh (7.41) tetapi oleh 

YP(1) #YP(JP 4-1) —0 
yang dapat ditulis ulang menjadi 


Set, H)-.x() —D(i) untuk i — JPt2 

rs 
dimana 

Cli,j)-0 untuk 2 jX JP danj- JP42 

Cli,1) — Ci, JP--1) —1 

D(i)-0 
Ini berarti bahwa untuk aerofoil tumpul, sistem persamaan setimbang yang harus diselesaikan 
adalah 


JP42 


» Ci i).x(i) —Dli) untuk 1 Sis JP42 (7.45) 
j-l 


dimana untuk 1Xi « JP 2 
Cli,1) — CPli,1) 
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Cli,j) - cP(i,j) -cO(i,j-1) untuk 25 j «JP 
Cli, JP --1) — COli, JP) 


wd 


(4) —-yP(j) untuk 1S js JP-41 
XJP 42)  Waero 

Di) - —2x/(XP(i).sina — YP/i).cosa) 
dan untuk i-JP--2 

Cj)-0 untuk2£ jS JP danj- JP-42 

Cli,1) — Cfi, JP 4-1) —1 

D(i)-0 
Sekarang marilah kita pelajari masalah metoda panel vorteks derajat tinggi untuk kasus 
kondisi batas Neumann. Komponen kecepatan arah-x, V,, dan arah-y, V,, yang terimbas di 


titik T karena keberadaan panel vorteks PO(j) dengan distribusi kekuatan linear diberikan 
oleh rumus-rumus berikut 


Vi) - CL) PP) C2). (rei) PG) (9.46) 
Vj) (CL 1.) Gi) C2 Ui). 1) PP) 9.41) 
dimana 





CL (T,j) -—A(T, j).cosB(j) ul La | Ja B(j) 








Jkoata A(T,j).sin B(j) 
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C2,(T,j) aa San RD. Ato eta 





poli) PT(T,j))) POlj 
ale PRO) al CM) TRTI) tp.) Lsingl: 
6 P0)" en) Poli ni PG) 





Na PR) | gal 
C2, (T,j) C PO) ml 


(at la 
Poli) ' 





Arah normal pada panel PO(i) adalah ni) dimana 


n(i) — sin B(j).x4-cosP(j).y 

cosB(j) - AXOP(i)/PO(i) 

sin B(j) - AYOP(i) /PO(i) 
Bila titik T dipilih berada pada tengah-tengah panel PO(i) maka komponen kecepatan normal 
pada panel PO(i) di titik T adalah 


Valisj) - V(i,i).cosB(j)— V,(i,i).sin BU) (7.48) 
Pemasukan persamaan (7.46) dan (7.47) ke (7.48) akhirnya memberikan hasil berikut 


Valisi) ——|(-At.co5(pti) —p1)) H B1.sin(gCi) —B(1)))-1PM5) 


#(-A2.co5(gli) -B(1)) -B2sin (BC) — PI))-oreli 1) TJ) 
(7.49) 


127 


dimana 
orki,j) 
al Meh) 


Bi-Ai,j) 





.NT,j) 





Aga PRT) i 3) TR/T,i) 
poli)” ter(r,i)) Poli) 





B2 - IRIT) ml OTT,i) | PRT) Ary) 
polj)) " leriT,i)) oli) 
Persamaan (7.49) dapat diolah lebih lanjut menjadi 


Vai) ---|((-A1 1 A2).co5(pMi) — PC) #(B1-4-B2).sin(B(G) -B(5)))- 7266) 


#(-A2.co5(Bli) - BC) B2.sin (BG) -(1)))-ve(i 1) 


V, (isi) — srllas.cosapti,i) — B3.sin AB(i,j)).vP() TE 
-(A2.cos AB(i,j) — B2.sin AB(i,j)). YP (jt 1)| 


dimana 





AT, j) 


2 MT) | OTT) | TR(T,j) 
sikat Tu (Ta) Poli) 


“MT, j) 





B3 RT) a Oa) aan 
poli)“ (PTIT,i)) Poli) 


AB(i, j) — B(i) —B(5) 


Syarat batas Neumann untuk kasus ini adalah 


sina.cosB(i) — cosa.sin PG) an Elon) Rt) cotiidanti #1)-0 


—1 
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Persamaan di atas dapat ditulis ulang menjadi 
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22) YP(j))—D(i) untuk 1S JP (7.51) 
dimana untuk 1 Xi « JP 
Cli,1) — CP(i,1) 
Cli,j) — CP(i,j) -CO(i,j-1) untuk 2 j« JP 
Cli, JP 4-1) —-COli, JP) 
( 


D(i) — sin x.cos Bli) - cosa.sin B4) - sin(P(i) — a) 
CP(i,j) — A3.cos AP(i, j) —B3.sin AFI, j) 


CO(i,j) - A2.cos AFli, j) -B2.sin AP(i, j) 
Karena T(i) adalah titik tengah panel PO(i) koordinatnya adalah 


| untuk 1Ki « JP 


xi) 1 (xe(i) « xp(i 1) 

Ti) - W (Ye(i) - ve(i--1)) 
Kondisi Kutta yang harus dipenuhi adalah 

YP(1) -YP(JP 4-1) —0 (7.52) 


Untuk kasus aerofoil berekor tumpul titik P(1) tidak berhimpit dengan titik P(JP-H1), jadi 
persamaan (7.52) dapat ditulis menjadi 


YP(JP 41) — —yP(1) (7.53) 
Pemasukan persamaan (7.53) ke (7.51) memberikan sistem persamaan berikut untuk 
Isis JP 


| untuk 1Xi« JP 


Sci) (57). P4) — Di) (7.54) 
ja 


dimana 


Cli,1) — CP(i,1) - coli, JP) 
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Cli,i) — CP(i,j) -C@(i,j 1) untuk 2 j «JP 


D(i) — 2x.sin(P/i) 5 a) 
Untuk kasus aerofoil berekor lancip, kondisi Kuta dapat dinyatakan sebagai persyaratan 
bahwa bisektor sudut ekor adalah garis arus atau kecepatan normal di titik B pada bisektor 
sudut ekor harus bernilai nol. Titik B dipilih berada sedikit di belakang titik ekor, yaitu 
koordinatnya adalah (XB, YB) dimana 


XB — XP(JP 4-1) --0.001 
YB -— tan (pi) 4 par) )o.oo1vetar 41) 


Titik B dapat ditentukan sebagai titik kontrol ke JP-#1, yaitu sebagai titik T(JP#1), dan 
persamaan (7.51) dapat diterapkan di (JP#1) titik konrol sehingga menjadi sebuah sistem 
persamaan yang setimbang. Data tambahan yang dibutuhkan adalah bahwa 

XT(JP 41) — XB 


YT(JP 4-1) — YB 


BIP 4-1) — V4 (B(1) -B(JP)) 

Sekarang marilah kita pelajari penerapan metoda panel untuk menganalisa aerofoil multi- 
komponen. Pada dasarnya ini tidak berbeda dari analisis aerofoil komponen tunggal, kecuali 
bahwa sistem indeks yang digunakan menjadi sedikit lebih rumit. Pertama-tama ditentukan 
bahwa aerofoil yang ingin dianalisa terdiri dari N komponen. Titik-titik pada permukaan 
aerofoil diberi indeks j yang nilainya meningkat secara berurutan. Titik ekor bawah 
komponen aerofoil ke n diberi indeks j-JEBIn) dan nilai j meningkat sesuai arah putaran 
jarum jam (lihat gambar (5.32) dan (7.32) sehingga akhirnya di titik ekor atas nilai indeks j 
menjadi JEA(n). Komponen aerofoil ke n mempunyai JP(n) panel, dimana 


JP(n) — JEA(n)- JEB(n) untuk 1 Kn N (9.55) 
Nilai-nilai JEB(n) dan JEA(n) dihubungkan oleh rumus-rumus berikut 
JEB(1) —1 | 


7.56 
JEB(n) — JEA(n—1)41 untuk2 LN en 
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Komponen aerofoil ke n diwakili oleh JT(n) titik-titik permukaan, dimana 


JT(n) — JP(n) 1 

JT(n) — JEA(n) — JEB(n)--1 (7.57) 

JT(n) - JEB(n #1) - JEB(n) 
Pada umumnya garis yang menghubungkan titik Pj) dan P(jt1) adalah panel PO). Ini 
berlaku untuk semua nilai j kecuali untuk j-JEA(n), 1X mX N-1, karena garis yang 
menghubungkan titik ekor atas komponen ke n dengan titik ekor bawah komponen ke (nt1) 
itu bukan suatu panel permukaan aerofoil (lihat gambar (7.3)). Metoda panel vorteks 
distribusi kekuatan seragam dengan kondisi batas Dirichlet dalam analisis aerofoil multi 
komponen dapat dijelaskan sebagai berikut. Variabel anu yang tak-diketahui dan harus 


dihitung adalah kekuatan vorteks y( j) di setiap panel untuk seluruh jumlah panel, JSP, 


dimana 


JSP - 5 y(n) (9.58) 


n-l 


Disamping itu kita juga tidak mengetahui nilai fungsi arus untuk masing-masing komponen 
aerofoil, yaitu @,ev, (1) dimana 1£ n £ N. Ini berarti bahwa secara keseluruhan ada JSP4N 


variabel anu yang tak-diketahui. Syarat batas Dirichlet dapat diterapkan di JP(n) titik pada 
permukaan komponen ke n, yaitu dari titik T(JEB(n)) sampai dengan titik T(JEA(n)-1), 
untuk masing-masing nilai n dimana 1X nXN. Syarat batas Dirichlet untuk komponen 
aerofoil ke n dapat ditulis sebagai berikut 


S( Sein) - Di) (9.59) 


k51 Vi-JEB(k) 
untuk JEB(n) «i « JEA(n) 1 
dimana nilai C(i,j)) dan D(i) adalah 
Clij) - p@li) -@R(i,i)in Oli, i) -PR(i,j) mPT(i,j) - TR(i,j)AG,j) 


Di) -— 2x XT) sina— YTi) cosa) 
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Bila dirumuskan bahwa 
Cli,JEA(k))-0 — untuk IKKEN dan k#n 
Cli, JEA (n)) 5-21 untuk k-n 
Y(JEA(k)) — 6... (K) untuk 1 KN 

maka persamaan (7.54) dapat ditulis menjadi 


JEA(n) 
5. cli,i)y(i) — bli) untuk JEBn) is JEA(n)-1 (7.60) 
j-l 
Kondisi Kutta untuk komponen aerofoil ke n adalah 
Y(JEB(n))-#y(JEA(n) —1) —0 
yang dapat ditulis ulang menjadi bentuk 


Hat, 1G) -—n(i) untuk ix JEA(n) (9.61) 


dimana 
C(i,j)-0 — untuk ISS JEA(N) kecuali j- JEB(n) dan j- JEA(n)-1 
Cli,j)-1 — untuk j-JEB(n) dan j- JEA(n)-1 
D(i)-0 
Sistem persamaan (7.60) dan (7.61) dapat digabung menjadi bentuk berikut. Untuk masing- 
masing nilai n, dimaan 1X n X N 
JeA(n) 


2 Cli, »)xG4) —D(i) untuk JEBn) sis JEA(n) (9.62) 


dimana untuk JEB(n) Kis JEA(n) —1 
D(i) - 2x(XT(i)sina — YT(i) cosa) 


Cli,i) - PO())-- @R(i, j)in Oi, j) - PR(i, j)InPT(i, j) - TR (i,j)A(,j) 
untuk semua nilai j, 1X jX JEA(N), kecuali j- JEA(k), 1ISkKEN 
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Cli,j)-0 untuk j-JEA(k), IXKXN kecuali k-n 

cli,j) -—21 untuk j- JEA(n) 
dan untuk i - JEA(n) (kondisi Kutta) 

D(i)-0 

Cli,j) -0 untuk semua nilai j, 1 jX JEA(N), kecuali j-JEB(n) danj — JEA(n)-1 
Cli,i) 1 untuk j- JEB(n) dan j- JEA(n)-1 
Perhatikan pula bahwa 

x(j) -y(i) untuk TK JEA(N), kecuali j-JEA(k), 1KkKEN 

x()) — @xera(K) untuk j- JEA(k), 1KKEN 
Sistem persamaan (7.62) terdiri dari JEA(N) persamaan dan melibatkan JEA(n) variabel anu 
sehingga merupakan suatu sistem persamaan setimbang yang dapat diselesaikan secara 
simultan dengan Eliminasi Gauss, L-U Dekomposisi atau metoda iteratif seperti Gauss- 
Seidel dan lain sebagainya. 

Metoda ini berlaku baik untuk komponen aerofoil berekor lancip ataupun tumpul dan 
titik kontrol atau titik kolokasi, T(i), boleh dipilih berada diujung panel, yaitu 


XT(i) — xP(i) 
untuk semua i, IX ix JEA(N) kecuali i- JEA(n), 1S n« N 
YT(i) — vP(i) 


(9.63) 


atau ditengah-tengah panel, yaitu 
xT(i) Hx16) #XT(i-1)) 


untuk semuai, 1SiS JEA(N) kecuali i - JEA(n), 1S n£ N 
yt(i) (XT) yr 1) 


(7.64) 
Sekarang marilah kita bahas kasus metoda panel vorteks distribusi kekuatan linier dengan 
kondisi batas Dirichlet. Untuk kasus ini titik kolokasi T(i) dipilih berhimpit dengan titik 
ujung panel P(i) untuk semua nilai i, dimana 1 XiX JEA(N), kecuali untuk komponen 
aerofoil "berekor lancip" dimana titik P(JEB(n)) berhimpit dengan titik P(JEA(n)), yaitu 
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XT(JEA (n)) - XP(JEA (n—1)) -0.99(xp/(JEA (n)) — XP(3EA (n—1))) 
(9.65) 

YT(JEA (n)) - YP(JEA (n—1))--0.99( YP(JEA (n)) — YP(JEA (n—1))) 
Syarat batas dirichlet untuk komponen aerofoil ke n, dimana 1 X n £ N, dapat ditulis sebagai 
berikut 


JEA(N)#N 
S cli,i)y(i) — Dj) untuk JEB(1) Sis JEAIN) (9.66) 
j- 
Rumus untuk C(i,j) dan Dff) adalah sebagai berikut 
untuk JEB(n) Ki JEA(n) dan 1K 1 N 
Di) —27.(XT(i).sina — YT(i).cosa) 
Untuk JEB(n) six JEA(n) dan 1£n« N 
untuk j — JEB(n) 
Cti,j) - CLi,j) - C2Gj) 
untuk JEB(n) #1sjs JEA(n) -1 
Ci) - CMi,j) - C2Gj) # C2G,j-1) 
untuk j — JEA(n) 
Cti,j) - C2ti,j-1) 
untuk JEB(N) 1 « j « JEA(N)#N 


Cli, JEA(N) 4-k) -0 untuk 1 KN bila kn 
Oli, JEA(N) Kk) --27 untuk 1 KN bila k-n 
Dalam rumus di atas C1(i,j) dan C2(i,j) adalah 
CKi, i) — PO(j) “OR, j).In OT(i,j) PR (i,j).InPT(i,j) - TR (i,j). Ali, j) 
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1 PO(i)-2PR(i,j)) PO) —PR(i,j) #TR(i,j) 





czi,j : an -nOT(i,j) 
PR ia) TR prt, j) PROd)-TR GI Ay) 





2PO/j) poli) 
berlaku untuk 1 Si « JEA(N) dan 1S j« JEA(N) kecuali j - JEA(N) dimana 1 1 XN . 
Kondisi Kutta untuk masing-masing komponen aerofoil dapat diberikan sebagai berikut 
Y(JEB(n)) #y(JEA(n))-0 untuk IK NN 
yang dapat ditulis ulang menjadi bentuk 
JEA(N)AN 
» Clii)(i) -D(j) untuk i — JEAMn) dan 1S SN 
ja 
dimana C(i,j) dan D(i) adalah 
D(i) —0 
C(i,j) - 0 untuk semua j kecuali j — JEB(n) danj - JEA(n) 
C(i,j)) -1 untuk j — JEB(n) dan JEA(n) 
Syarat batas Dirichlet (7.66) dan kondisi Kutta (7.67) dapat digabung menjadi satu sistem 
persamaan setimbang, terdiri dari JEA(N)#N persamaan dan melibatkan JEA(N)#N variabel 


anu, sebagai berikut 
Untuk 1S i« JEA(N)#-N 


JEA(N)#N 


2 Clii)(i) - Di) (9.68) 


dimana C(i,j) dan D(i) diberikan sebagai berikut 
Untuk 1S ns N 
Untuk JEB(n) Si JEA(n) 


D(i) 2x. (XT(i).sina — YT(i).coso) 
Untuk 1S ms N 
Cli, JEB(m)) - Cili, JEB(m)) - C2(i, JEB(m)) 
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YT,»(i) 5 XT(i).tana — » C,(ij).1G) untuk 1S is JTA (8.28) 


ja 
Pengurangan persamaan (8.27) oleh (8.26) dan setelah dilakukan pengolahan lebih lanjut 
akhirnya dapat memberikan hasil berikut 


VI YA Gali) li, »))() untuk KIS ITA (8.29) 


Persamaan di atas dapat juga ditulis ulang menjadi bentuk 
YToa(i) — YT, (i) #AYT, (i) untuk 1 Si JTA (8.30) 


dimana AYT, (i) adalah suku perbaikan atau koreksi, yaitu 


AYT, (i) Hai) —C,(i,j))-y(i) untuk 1 ix JTA (8.31) 


Persamaan atur untuk aliran potensial tak-mampat adalah persamaan Laplace, yaitu suatu 
persamaan diferensial parsial derajat 2 yang bersifat eliptik dan linear. Persamaan eliptik 
dapat diselesaikan secara numerik dengan mengunakan matode taksiran beda berhingga 
(TBH) atau finite difference approximation. Dalam pendekatan numerik dengan metode 
TBH ini solusi diperoleh secara iteratif dengan menggunakan rumus iterasi yang berbentuk 
sama dengan rumus (8.29). Untuk mempercepat laju konvergensi, rumus iterasi tersebut 
diperbaiki dengan menerapkan konsep SOR atau Successive Over Relaxation, yaitu 
perbaikan beruntun berlebihan. Kuantitas AYT, (i) dalam rumus (8.29) adalah koreksi atau 
perbaikan nilai YT(i) dari iterasi satu ke iterasi berikutnya. Walaupun konsep SOR 
sesungguhnya berlaku untuk metoda TBH yang diterapkan pada -persamaan Laplace, tetapi 
diperkirakan bahwa konsep tersebut berlaku juga untuk metoda panel yang diterapkan pada 
persamaan Laplace. Dengan pertimbangan tersebut maka rumus iterasi (8.29) dapat dicoba 
diperbaiki sebagai berikut 
XT (1) 5 YT, (i) to. AYT, (i) (8.32) 

dimana @ adalah faktor SOR yang untuk metoda TBH bernilai di jajaran berikut 


1c0x2 
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Pada awal bab ini dinyatakan bahwa data masukan untuk masalah perancangan adalah 
distribusi koefisien tekanan, Cp disepanjang permukaan aerofoil. Dalam prakteknya 
permukaan aerofoil diberikan oleh koordinat titik-titik P yang berada pada permukaan 
aerofoil, yaitu (XP())YPG)) untuk 1X jX JTA. Pada umumnya nilai indeks j ditentukan 


meningkat sesuai arah putaran jarum jam dan titik P(1) adalah titik ekor permukaan bawah 
sedangkan P(JTA) adalah titik ekor permukaan atas aerofoil (lihat gambar (8.1)). Jumlah 
titik aerofoil, JTA, adalah bilangan ganjil. Ini berarti bahwa indeks untuk titik titik hidung 
atau titik aerofoil di ujung depan, bernilai 

JH — Y (JTA #1) (8.33) 
Data bentuk aerofoil diberikan sedemikian rupa sehingga 

XP(j) — XP(JTA —j--1) (8.34) 
Distribusi Cp yang diberikan sebagai data masukan dalam prakteknya diberikan sebagai 
sejumlah pasangan data (XP(j), Cp(j)) untuk 1S jX JTA. Distribusi kekuatan vorteks, 


y(i), kemudian dapat dihitung dengan rumus (8.2) dan (8.3) yaitu 


1) -J1-Cp(j) (8.35) 
dimana y(j) bernilai positif untuk permukaan bawah sedangkan untuk permukaan atas 
aerofoil y(j) bernilai negatif. Permukaan bawah dan atas adalah permukaan-permukaan 
yang dibatasi oleh ekor dan titik stagnasi depan dimana y — 0. Titik stagnasi depan tersebut 


berada di panel PO(JS) dimana JS tidak boleh bernilai lebih daripada JH. Untuk aerofoil 
simetris pada sudut serang « — 0, titik stagnasi depan berhimpit dengan titik hidung atau 


y(JH)-— 0, yaitu Cp(JH) — 0. Untuk aerofoil lengkung atau aerofoil simetris pada sudut 
serang a » 0, titik stagnasi depan berada di panel belahan bawah di belakang titik hidung. 


Sebelum tanda untuk y(j) dalam rumus (8.34) dapat ditentukan kita harus menentukan 
panel JS dimana titik stagnasi depan berada. Nilai mutlak y(JS) dan y(JS 1), sesuai 
penalaran wajar, adalah 2 nilai mutlak y(4) yang paling kecil. Pertama-tama kita harus 


memeriksa nilai indeks j dimana y( j) bernilai minimal. Prosedur yang perlu dilakukan adalah 


sebagai berikut. Pada awalnya ditentukan bahwa 
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Yin SY (JH) 
Suatu prosedur saringan atau sortir, kemudian dilakukan sebagai berikut 
Untuk j — JH-1, JH-2, JH-3, . . .dan seterusnya 
Bila Yi 21 j) maka nilai y,,,, diperbarui, yaitu 

Ymin SU) 
Persyaratan di atas akan selalu dipenuhi sampai akhirnya y,,i, adalah benar-benar nilai 
minimal dari y(j). Pada saat pertama kali persyaratan y,,,, » y(j) tidak terpenuhi, nilai j 
adalah JS-1, jadi dicatat bahwa JS-j-1. Selanjutnya penalaran wajar akan menyatakan 


bahwa nilai y(j) terkecil kedua boleh berada di j-JS-1 atau j-JS#1. Untuk menentukan nilai 
j dimana y(j) adalah nilai y terkecil kedua, kita dapat melakukan pengujian berikut 
Bila Iy(3S —1)| Iy(3S 4-1)| maka y(JS— 1) adalah nilai y terkecil kedua. Untuk kasus ini 


titik stagnasi depan berada di panel PO(JS-1) dan y(i) bernilai positif hanya untuk 
ISjxJS-1 dan bernilai negatif untuk JSSJSXJTA. Sebaliknya bila 
(xtas- 1) » Ix(as--1) maka y(JS-4-1) adalah nilai y terkecil kedua. Untuk kasus ini y(j) 
bernilai positif untuk 1S jS JS dan bernilai negatif untuk JS#-1Xj« JTA. Dengan 
pengertian tersebut di atas, prosedur atau algoritma perancangan bentuk aerofoil dengan 


matoda panel dapat dijabarkan secara lebih rinci sebagai berikut. 
Pertama-tama distribusi (XP(j), Cp/(j)) untuk 1X jX JTA diberikan sebagai data 


masukan Nilai-nilai y(j) dengan tanda yang benar kemudian dihitung seperti dijelaskan di 
atas. Titik-titik kolokasi T(j) ditentukan berhimpit dengan ujung panel, jadi pada umumnya 
XT) - XP()) 
kecuali untuk j-JTA, dimana 
XT(JTA) — XP(JTA)—0.001(XP(JTA) - XP(JTA —1)) 
Kemudian nilai-nilai ordinat YP(j) ditebak, misalnya dengan cara menggunakan bentuk 
aerofoil NACA-0012 (atau bentuk aerofoil sembarang lainnya) sebagai bentuk aerofoil 
tebakan awal. Nilai-nilai YT(j) tebakan kemudian dihitung, yaitu 


YTH) - YP(H) 

YT(JTA) — YP(JTA) —0.001( YP(JTA) — YP(JTA —1)) 
Selanjutnya nilai-nilai koefisien Cfi,j) tebakan dihitung dengan rumus-rumus (8.6), (8.7), 
(8.8), (8.9) dan (8.10) dimana 
poci) - (rei) xD) (ve) — yet) )” 
POC)-PR(i,j) - (XP(j--1)-— xP()) (xi) — xP(1))-“( XP(3--1)-— yP(1)) xi) — yel) 
POC)-TR (i,i)  (XP(i--1)-— xp(1))( xar(i) — ve(1))— (xe(i--1)— xe(1))(xar(i) — xe(1)) 


POXi).OR(i, j) - PO(j).PR(i, j)- PO(j) 


Pri) - (xxi) - xe(3)) (v1) ve)” 


PA 


OT, ) - (XTG) — XP) Avi) — ve) 





AGsid-tan POD.TR GH) 
PT jY -PO(i)-PRG, j) 
Rumus-rumus diatas dan rumus (8.9) dan (8.10) berlaku untuk kasus dimana titik T(1) tidak 
berhimpit dengan titik P(j). Apabila T(i) berhimpit dengan P(j)) maka rumus (8.9) dan (8.10) 
disederhanakan menjadi bentuk berikut 
CK) - PO(i)-(1- In PO(i)) 
CZi,i) -0.sPOi).(0.5— InPOXi)) 
Apabila T(i) berhimpit dengan O(j) atau P(j-1) maka 
CK Lj) - PO(j).(1-In PO())) 
CAj-#L,j) — 0.5PO(j).(1.5-InPO()) 
Perbaikan awal nilai ordinat YP(j) kemudian dilakukan sebagai berikut. 


JTA 
Nilai awal » C(i,j).y(4) dihitung dan disimpan sebagai data F(i) berikut 
j- 
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Fi) - S Cd aH untuk 1 Si JTA 


ja 

Perbaikan awal ordinat YT(i) dan YP(j) dilakukan sebagai berikut 

YT(i) - XT(i).tana —F(i) untuk 1 Si JTA 

YP(j) - YT(j) untuk 1 jK JTA -1 

YP(JTA) — (YT(JTA) — 0.001. YP(JTA —1)/0.999 
Nilai YT(i) dan YP) selanjutnya dihitung secara iteratif sebagai berikut 
1. Hitung Cfi,j) dengan rumus (8.6), (8.7), (8.8), (8.9) dan (8.10) atau bentuk khususnya 

bila i-j dan bila i-j--1 

2. Nilai F(i) lama disimpan, misalnya sebagai AYTi) 
3. Nilai baru F(i) dihitung, yaitu 


Fi) — S Cadar untuk 1 ic JTA 
ja 


4. Nilai AYT(i) yang sebenarnya, yaitu AYT(1) - Fi)saru " F(i)uma Gihitung , jadi 
AYT(i) -F(i)-AYT() untuk 1S i£ JTA 
5. Uji konvergensi 
Bila JAYT(i)| » & untuk salah satu nilai, proses iterasi belum konvergen dan nilai YT(i) 
dan YP/(j) diperbaiki, yaitu 
YT(i) - YT(i) to. AYT(i) untuk 1 Li JTA 
YP(j) - YT(j) untuk 1S JK ITA-1 
YP(JTA) — (YT(JTA) — 0.001. YP(JTA —1)/0.999 
Iterasi kemudian diulang lagi mulai dari butir i. 
Bila (AYTG)| Se untuk semua nilai i, proses iterasi telah konvergen. Nilai YT(i) dan YPG) 


dihitung untuk terakhir kalinya sebagai berikut 
YT(i) - YT(i) #AYT(i) untuk LSI X JTA 
YP(j) - YT(j) untuk 1K j JTA-1 


YP(JTA) — (YT(JTA) — 0.001. YP(JTA —1)/0.999 
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Iterasi dihentikan karena ordinat YP(j) yang sesungguhnya sudah diperoleh dengan derajat 
keakuratan sebesar £ . 

Metoda panel Vortex dengan syarat batas Neuman juga dapat diterapkan untuk merancang 
bentuk aerofoil. Kalau dalam syarat batas Dirichlet permukaan aerofoil adalah sebuah garis 
arus tertentu maka dalam syarat batas Neuman komponen kecepatan arah normal pada 
permukaan aerofoil harus bernilai nol. Permukaan aerofoil diwakili oleh panel-panel dan oleh 
karena itu titik kolokasi untuk kasus ini sebaiknya dipilih berada ditengah panel. 

Dari bab 7 diketahui bahwa komponen kecepatan arah-x atau arah-y, yaitu u dan v, dititik 
kolokasi T(i) diberikan oleh rumus-rumus berikut 


JTA-1 


an MP en 
uli) — cosa aa LAGDAD 


, 1 ITA PRG,j)) , OT(i,j)) PR,j) i —1G 
at (w POC) “"PTG) POM agro Na 











PA Ta n LG) 
v(i) — sin ja DI PT) yG) 








La RB), TG) TRI) 
Tan T PO) PT) POM) Gp Jna- 1) 


Rumus untuk komponen kecepatan normal adalah 
Vnfi) ——u(i).sin B(i) -v(i).cos B(i) 
dimana 
sin B(i) — AYOP(i)/PO(i) — (XP(i #1)— YP(i))/POi) 
cos B(i) - AXOP(i)/PO(i) — (XP(i-#1) — XP(i))/PO(i) 
Rumus-rumus di atas setelah diolah secara aljabar yang cukup rumit dan berkepanjangan 


akhirnya dapat ditulis menjadi bentuk berikut 


JTA-1 


Vnfi) — —sin PG).co.#.4 cosPG)-sina 4-7 HC, 1) 261,1) GD 1G) 


ja 
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dimana 


CX, j) — Ali, j)sin PA cos) (8.36) 


C2,j) ca cosB(i) 2 PR(i, j).cos B(i) — TRGi, j).sin B(i) In OTG, ) 








PO(j) PTKi,j) (837) 
PR, H-sinB()— TR, H-cosB(i) AGi,i) | 
PO) 
Rumus Vnfi) di atas dapat juga ditulis sebagai 
Vnfi) 4-sin B(i).cosa — cos B(i). sin Ya (8.38) 
ja 
dimana 
C(i, 1) — CL, 1) - C2, 1) 
Cfi,j) — CKi, j) -C2(i,j)-C2(i,j 1) untuk 2 jX JTA-1 (8.39) 
C(i, JTA) —-C2(i, JTA) 
Khusus untuk kasus i-j, XT(j) — 0.5(XP/(j)-- XP(j -1)) 1x j JTA 
CXj,j) —2n.sin P(j) 
40 
ne apain LAN is 


Syarat batas Neuman, Vn — 0, dapat ditulis sebagai 


sing(i).cosa — c0sB(.sina4--Y CG, )-1) untuk 1SiS JTA-1 (8.41) 
ja 


Persamaan untuk i — JTA adalah kondisi Kutta, yaitu 

1) #Y(JTA) —0 (8.42) 
yang harus dipenuhi oleh distribusi y(j) yang diberikan. 
Mengingat bahwa 

sin B(i) - AYOP(i)/PO(i) 

cos Bi) - AXOP(i)/PO(i) 
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maka syarat batas Neuman dapat ditulis ulang menjadi 


AYOP(i) — AXOP(i). tang, ta ro $ cd, D-10) (8.43) 


aj 51 


Koordinat titik-titik permukaan aerofoil diberikan oleh rumus berikut 


YP(1) —K (suatu konstanta) | (8.11) 


YP(i #1) — YP(i) #AYOP(i) untuk 1 Li JTA -1 
Bentuk aerofoil diberikan oleh distribusi nila AYOP(i) sedangkan nilai K hanya 
menentukan letaknya di sepanjang sumbu-y. Pengamatan ini menunjukkan bahwa dalam 
metoda panel dengan kondisi batas Neuman rumus (8.42) adalah rumus yang berfungsi 
seperti rumus (8.23) dalam metoda panel dengan kondisi batas Dirichlet, yaitu dapat 
digunakan untuk memperbaiki tebakan bentuk aerofoil. 


Pada awalnya bentuk aerofoil atau kordinat (xp(i), xp, (1)) ditebak. Koordinat titik 
kolokasi, T(i), kemudian dihitung sebagai berikut untuk 1 i X JTA -1. 


XTKi) — X (xeli) #xP(i--1)) 


YT,(i) — V4 (YP,(i) 4 Ye,(i--1)) 
Perhatikan bahwa nilai-nilai XP(i) dan XT(i) diberikan sebagai data dan tidak perlu 
diiterasikan. Nilai YP(i) juga diberikan, misalnya YP(i)-0, jadi tidak perlu ditebak ulang atau 
diiterasikan. 
Dengan pengertian tersebut maka hanya nilai-nilai ordinat YP(##1) dan YT(i) untuk 
1Sis JTA-1 saja yang perlu diiterasikan. Rumus-rumus untuk YP(i-1) dan YT) dapat 
ditulis ulang menjadi 


YP(i #1) — YP(i) -AYOP(i) 

YT(i) — XP(i)--0.5. AYOP(i) 
Setelah nilai YP,(i-1) dan YT,(i) ditebak, kemudian nilai koefisien pengaruh C, (i,j) 
dihitung menggunakan rumus-rumus (8.35), (8.36), (8.38) dan (8.39). Rumus (8.42) 
kemudian diterapkan untuk menghitung AYOP, (i) , yaitu 


| ntu ISIS JTA-1 
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JTA-1 


AYOP, (i) — AXOPMi). Magna Pat 3 @).7G) (8.45) 
Nilai-nilai YP(i#1) dan YT(i) kemudian diperbaiki Th berikut 


YP, (2) — AYOP, (1) 
YP, (1-1) — YP, (i)--AYOP, (i) untuk 2 Si JTA-1 


YT, (1) —0.5.AYOP, (1) 

YT, (i) — YT, (i—1)-#0.5.AYOP, (i—1) untuk 2S JTA-1 
Nilai C,(i,j) kemudian dihitung sebagai fungsi YP, (j) dan YT, (1) dan diterapkan untuk 
menghitung AYOP, (i) yaitu 


F9 kar 


AYOP, (i) - AXOP(i).tana 4d Ng @&).G) 


Rumus umum untuk iterasi ke-n adalah 


AYOP, (i) — AXOP(). tan Om Oa Co 3) 


cosa. 15 
yang dapat juga ditulis sebagai berikut 
AYOP, ,, (i) — AYOP, (i) #5, Y(i) (8.46) 
dimana 
P JTA-1 P JTA-1 
2G AU PLAN aa NADA 


2x 
Konsep SOR dapat Pa untuk memperbaiki rumus (8.45) menjadi 
AYOP, ,, (i) — AYOP, (i)t@.5,. Y(i) (8.47) 
Algoritma perancangan bentuk aerofoil berdasarkan metoda panel Vortex distribusi linier 
dengan syarat batas Neuman dapat dijelaskan sebagai berikut. 


A. Distribusi koefisien tekanan diberikan sebagai nilai-nilai (XP(j)))CPG)) dan 
dibaca sebagai data masukan untuk is j«s JTA. 


Nilai-nilai y (j) dihitung dengan rumus 


164 


1) - -J1- Cpli)” 


Nilai-nilai XTG) dihitung dengan rumus 
XT(j) —0.5(XP(j) - XP(j--1)) 
Nilai YP(1) ditentukan, yaitu YP(1)-0 
Bentuk aerofoil ditebak. Nilai tebakan awal YP(j) ditentukan dan YT(j) dihitung 
dengan rumus 
YT() —0.5(XP(j)  VP(j--1)) 
Nilai koefisien pengaruh C(i,j) dihitung. 
Nilai tebakan awal AYOP(i) dihitung dengan rumus 
AYOP(i) — (XP(i--1)— XP/i)).tana 4 F(i) 
dimana F(i) telah dihitung sebelumnya, yaitu 
P i JTA-1 Ta 3 
TO CG 


F@)— 
@ 21.c0osa 5 
Perbaikan awal bentuk aerofoil dilakukan. 


Nilai-nilai YP(#1) dan YT(i) dihitung dengan rumus 
YP(i-41) — YP(i) #AYOP(i) 
YT(i) — YP(i)--0.5. AYOP(i) 
Prosedur iterasi untuk memperbaiki bentuk aerofoil dimulai dan diteruskan 
sampai konvergen. 
Nilai C(i,j)) dihitung sebagai fungsi nilai terakhir koordinat (XP),YPG)), 
(XTM),YT@)). 
Nilai terakhir F(i) yang disimpan, yaitu 
SY) -F@) 
Nilai F(i) yang baru dihitung 


PO) Wi 


FO ena ACA 


Nilai SY(i) yang sebenarnya dihitung, yaitu 
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sXli) - Fi) -SY(i) 
Galat iterasi dihitung, yaitu 
galat(i) — (SY (i)| 
Galat maksimal dicari 
galat maksimal — nilai galat (i) terbesar 
Uji konvergensi 
Bila galat maksimal »& , iterasi belum konvergen 
Nilai AYOP(i) terbaru dihitung 
AYOP(i) —- AYOP(i) -o.5Y(i) 
Nilai YP dan YT terbaru dihitung 
YP(i-1) — YP(i) #AYOP(i) 
YT(i) — YP(i) t0.5.AYOP(i) 
Kembali ke awal iterasi 
Bila galat maksimal « &, iterasi sudah konvergen 
AYOP(i) — AYOP(i) tSY(i) 
YP(i--1) — YP(i) #AYOP(i) 
YT(i) — YP(i) 0.5. AYOP(i) 


Selesai 
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untuk JEB(m) 41 Si « JEA(m)-1 
Cli,)) — cx(i,j) - c2(i,j)-c2(i,j —1) 
Cli, JEA (m)) - C2(i, JEA (m) —1) 
untuk 1S MX N 

Cli, JEA(N) --m) 0 bila mn 


Cli, JEA(N)--m) ——217 bila m-n 
untuk 1S mn KN 
i-JEAIN) tn 
Di) -0 
untuk 1 £ j « JEA(N) 4 N kecuali j-JEB(n) danj — JEA(n) 
Cfij) —0 
C(i,JEB(n)) — 1 
Cti,JEA(n)) — 1 
Perlu diingat bahwa dalam rumus (7.68) 


x(j) -y() untuk 1 « js JEAIN) 

x(j) — Yaer (1) untuk j- JEA(N) #n dan 11X N 
Penerapan metoda panel vorteks distribusi kekuatan linear dengan kondisi batas Neumann 
dapat dijelaskan sebagai berikut 
Untuk kasus ini titik-titik kolokasi, T(i), boleh dipilih berada di ujung atau di tengah-tengah 
panel PO(i). Apabila T(i) dipilih berada di tengah panel PO(i) maka arah normal diberikan 
oleh rumus berikut 


PN AN 


n(i) — —sinBli).xt-cosBli).y (9.69) 


cosBli) — Na (9.70) 


sinB(i) - AYOP(i)/PO(i) 
Perlu diperhatikan bahwa P(i) boleh bernilai positif atau negatif dan nilai mutlaknya berada 


di selang 


136 


Os Bi) xx 7.71) 


Untuk kasus dimana T(i) berhimpit dengan P(i), arah normal diberikan oleh rumus berikut 


nli) — —sinOfi).xx-cos0li).y 0.72) 
O(i) — (1—1(i)).BGi-1)--1).P(1) Maa 

r(i) - poli-1)/(roli-1)-ro(i)) | 
Rumus (7.73) berlaku untuk semua titik permukaan komponen aerofoil kecuali untuk titik 


ekor atau i — JEB(n). Untuk titik ekor, titik kolokasi T(i) dipilih berada di tengah-tengah 
titik ekor bawah dan atas, yaitu 


XT(JEB(n)) Hxe(sentu)) #XP(JEA (n))) 
(9.74) 


YT(JEB(n)) - Hxe(sen(n)) #YP(JEA(n))) 


dan sudut @ diberikan oleh rumus berikut 


H(p(oeata)-plen(a)-x) 079 


Perhatikan bahwa untuk kasus metoda ini, prosedur penyelesaian masalah untuk aerofoil 


O(JEB(n)) — 


berekor lancip dan tumpul itu sama saja. Prosedur dengan titik kolokasi, T(i), berada di 
tengah-tengah panel itu jauh lebih sederhana dibandingkan dengan untuk T(i) berhimpit 
dengan ujung panel P(i), tetapi ada kalanya kita harus menggunakan kondisi T(i) berhimpit 
dengan P(i). Persamaan atur untuk masalah ini dapat dijelaskan sebagai berikut 
Untuk setiap komponen aerofoil atau masing-masing nilai n dimana 1X n X N 

Untuk JEB(n) sis JEA(n)-1 (syarat batas Neumann) 


Seth) -DI)) (4.75) 


jA 
dimana 


Di) - 2x.sin (Oli) -a) 
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untuk j — JEBIn) 
Cli,j) — Cili,j).cos(eli) — P(j)) -— czli,i).sin(O(i) —P(1)) 
untuk JEB(n) #1£ jx JEA(n)-1 


cli,j) - Cili,i).cos(@(i) — B(j))-— cz(i,i).sin(O(i) —B())) 
#Cafi,j—1).sin(O(i) — B(j —1))-C3(i,j —1).cos(0(i) — Bli —1)) 
untuk j — JEA(n) 

Cli,j) - Cafi,j —1).sin(O(i) — B(j —1))- C3(i,j —1).cos(O(i) — Bli —1)) 

Untuk i-JEA(n) (kondisi Kutta) 

JEA(N) 
» iii) — oli) (7.76) 
j-l 

D(i) -0 

untuk j — JEBIn) 

C(,j) — 1 

untukJEB(n) #1£ js JEA(n)-1 

C@,j) -0 

untuk j — JEA(n) 

C@,j) 1 

Rumus untuk CNi,j), C2(i,j), C3(i,j), dan C4(i,j) adalah 


cili,j) — 





OR) (eri) TR)... 
og) Te) Po) AG 


c2ti,i) — 





TR(4j) | ( OT), ORG) 1. . 
Poli) aan) Po) Aa) 


Ai) 





C3(i,j) —1- 


PR(i,) wi Me) 
roti) (tertij)) poli) 
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Ai Toy “erGa) 7 reg) 

Sistem persamaan (7.75) dapat digabung dengan (7.76) membentuk suatu sistem 
persamaan setimbang yang melibatkan JEA(N) variabel anu dan terdiri dari JEA(N) 
persamaan. Sistem persamaan setimbang tersebut dapat diselesaikan secara simultan dengan 


TR(i,i) | Tea) at) Ati) 


salah satu metoda standar seperti Eliminasi Gauss dan lain sebagainya. Hasil yang diperoleh 
adalah nilai y di setiap titik P(j) untuk 1 « j X JEA(N). Apabila YP(j) adalah nilai y di 
titik P()) maka y di titik sembarang, S, dapat dihitung sebagai berikut 

SAR rei) say a77) 
dimana PS adalah jarak dari P(j)) ke S yang berada pada panel PO(j). Mengingat bahwa 
kecepatan aliran, V, di titik S itu sama dengan besaran y di titik S (lihat persamaan (7.21) 


dan (7.22)) maka koefisien tekanan di titik S adalah 
Cp(s) -1-—yS? (9.78) 
Koefisien tekanan rata-rata pada panel PO(j) dapat dihitung sebagai berikut 


Poli) 
Cpli)- | cpls).as7poli) 





polj) 2 
Cpti)- J (-Latstataon eat) lawreta 
Cpli) 1-3 (7) 4126) YP) rei) (279) 


Gaya yang beraksi pada panel PO(j) adalah gaya non-dimensional C, (j) dimana 


Crli) - Cpli).POi) (7.80) 
Gaya tersebut dapat dianggap beraksi di titik M(j) dimana jarak dari P(j)) ke M(j) adalah 
PM(G), yaitu 
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poli) 
PM(i) -| Jewel) fe (1) 


1 6-(netif 2np(i). reli 41) -3ye(i si) 
2 6 Ali) #ye(i) meli ta) 1yeli 1) 


Perhatikan bahwa seandainya YP(j) - YP(j-1)-y(j) persamaan (7.79) dan (7.81) 


PM(j) 








(7.81) 


disederhanakan menjadi 


Cpli) -1-—x(i) (9.82) 


PM(i) - -POfi) 283) 


yaitu sama dengan untuk kasus distribusi kekautan seragam. Koefisien gaya angkat dan gaya 
tahan dapat dihitung dengan rumus berikut 


C- -$ Cp(s).cosB(s)-ds/tali busur 


C,- -$ Cp(s). sinP(s). ds/tali busur 


atau 


JEA(n)-1 
- Yi) .cosPlj Ji busur 


i-JEB(n) 


Ca(n) - Bg sr (4)-sin B(j Jai busur 


i-JEB(n) 
dimana C,(n) dan C,(n) adalah koefisien gaya angkat dan gaya tahan untuk komponen 


aerofoil ke n Rumus-rumus di atas dapat diuraikan lebih lanjut menjadi bentuk berikut 


-1 


Eh Cp(j).AxoP(j) Jaa busur (7.84) 


i-JEB(n) 
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JEA(n)-1 


- X Cplj).Avorli) Jami busur (7.85) 


i-JEB(n) 


Koefisien gaya angkat dan gaya tahan untuk aerofoil secara keseluruhan adalah 


C-»XGln) (7.86) 


C, - » aa) (7.87) 


Gaya-gaya yang beraksi pada masing-masing komponen aerofoil juga menimbulkan momen- 
momen gaya relatif terhadap titik-titik engsel tertentu. Bila gaya resultante atau gaya total 
yang beraksi pada komponen ke n beraksi pada titik pusat gaya, PG(n),maka momen di 
PG(n) tentu saja bernilai nol. Ordinat PG(n) boleh ditentukan bernilai nol sedangkan 
absisanya dihitung sebagai berikut 

Momen di titik PG(n) adalah 


"SU On)-axor6) em) — rata) 


i-JEB(n) 


#Cpli).AYOP(j)(XPM(j) — xPG/(n))| (7.88) 


dimana (XPM(j), YPM()) adalah koordinat titik PM), sedangkan (XPG(n), YPG(n)) adalah 
koordinat titik PG(n). Koordinat PM) dapat dihitung dari rumus (7.81) yaitu 


PM(j) 
poli 
PM/(j 
Poli) 
Pemasukan (7.89) ke (7.88) dan mengingat bahwa YPG(n)-0 akhirnya dapat memberikan 
hasil berikut 


y: .AXOP(j) 
(9.89) 
) .AYOP/j) 


XPM(j) — XP(j) 








YPM(j) — YP(j)- 
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"S1 @IAYOR(j).xeM) - AxOP(j).vPM((j)| 
XPG (n) — SE) 





JEA(n)-1 


» Cpli).Avor(i) 


i-JEB(n) 


XPG/(n) - "S1 @Iaxor(i).xe(i) - AXOP(j). vP(j)| Ke 


i-JEB(n) 
(7.90) 
dimana tb adalah panjang tali busur aerofoil seutuhnya. Nilai momen di titik engsel EM 
dengan koordinat (XEM, YEM) diberikan oleh rumus berikut 


Momen di EM — P3 Cpli)favor(i).(xem(j) — xem) — Axor/(i).( vem(|) — vm) 


Persamaan di atas dapat diolah lebih lanjut menjadi 


Momen di EM - | C,(n).(XPG(n) — XEM) - YEM.C (j) |-tb (7.91) 
Koefisien momen di EM untuk elemen ke n adalah 
C,(EM,n) 2 Simone EM 
tb 
C,(n).(XPG(n) - XEM) -— YEM.C, (j 
C, (EM,n) - s(n) ( (n) ) (5) (7.92) 





tb 
Titik pusat tekanan aerodinamika untuk seluruh aerofoil adalah titik PA dengan koordinat 
(XPA,YPA) yang nilainya dapat dihitung sebagai berikut 


cun 5. |G()-axort) vem(i)— ra) ja 


#C, (i).AYOP(j).(XEM(j) — xpa))! -0 


n-l 
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Ordinat YPA boleh ditentukan sesuka hati dan untuk mudahnya dipilih bernilai nol sehingga 
koordinat PA adalah (XPA,0) dimaan nilai XPA adalah 


pa GP LAYER) Max) 








XPA n-1 | j-JEB(n) — 
2 jpn G6.avor)| 
atau 
XPA - Y'XPG/n).C,(n) S cln) (7.94) 


Karakteristik aerodinamika suatu aerofoil biasanya diberikan sebagai informasi tentang 
koefisien-koefisien aerodinamika, yaitu koefisien gaya angkat, koefisien gaya tahan dan 
koefisien momen. Nilai-nilai koefisien aerodinamika tersebut dapat dihitung menggunakan 
metoda panel dan rumus-rumus yang telah dijelaskan di atas 
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LAMPIRAN 7.1 


I. Metoda Panel Vorteks distribusi seragam dengan kondisi batas Dirichlet. 
I.1. Aerodinamika berekor lancip. 

Titik kontrol, T, dipilih berhimpit dengan titik ujung panel, P. 

Untuk 1xixs JP 

xTli) - xP(i) 

YT(i) — vP(i) 
1.2. Aerofoil berekor tumpul. 


Titik kontrol, T, dipilih berada di tengah-tangah panel PO. 
Untuk 1Xix JP 


XT) — U (xeki) - XP(i--1)) 
yT(i) - W (YP(i) - ve(i--1)) 
Persamaan Atur 


JP 


» Ci i).x(i) —Dli) untuk Sis JP #1 
ja 

dimana 

Untuk 1Sis JP (syarat batas Dirichlet) 
D(i) - 21(XT(i).sina — YT(i).coso) 
Untuk 1xixs JP 


Clij) - PO(j)-- OR, j).In OT(i,j) — PR (i, j).InPT(i,j) - TR(i, j). Ali, j) 
x(j) — y(i) 


Untuk i-JP-1 


Cli,j) - —2x 
X(4) - Waero 
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Untuk i-JP#-1 (kondisi Kutta) 


D(i)-0 

Untuk j — 1 danj — JP 

cli,j) 1 

Untuk 2S j X JP-—1 danj — JP 41 
cli,j)-0 


II. Metoda Panel Vorteks distribusi seragam dengan kondisi batas Neumann. 
Titik kontrol, T, dipilih berada di tengah-tengah panel PO. Bila semua titik tengah 
digunakan, akan memberikan sebuah sistem persamaan dengan ketentuan berlebihan 
atau "an over determined system of eguations". Oleh karena itu salah satu titik kontrol, 
di panel PO(JP), tidak digunakan, kasus ekor lancip dan tumpul sama saja 
Untuk 1S is JP-1 


xi) « W (Xi) — xP(i 1) 
vTti) - 1 (ve(i) - ve(i--1)) 


Persamaan Atur 


Cli,j).x()) —D(i) untuk1x ix JP 


Ms 


1 


dimana 
Untuk 1Six JP-1 (syarat batas Neumann) 


D(i) —21.sin(B(i)—a) 
Untuk 1£ ix JP 





Ci) - Ha Jodo -B(H))--Ali,i).sin(8G) -B6) 


Untuk i-JP (kondisi Kutta) 
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D(i) -0 
Untuk j — 1 danj — JP 
cli,j) -1 
Untuk 2 j £ JP-—1 
cfii) -0 
III. Metoda Panel Vorteks distribusi linear dengan kondisi batas Dirichlet. 
Titik kontrol, T, dipilih berhimpit dengan titik ujung panel, P. Untuk aerofoil lancip titik 
P(1) berhimpit dengan titik P(JP--1), sehingga untuk kasus ini y(JP 1) — y(1) atau 
14) (IP 41) —0 


Persamaan Atur 
JP-4-2 


HC,).xG) —D(i) untuk 1 Lis JP 42 
jl 


dimana 
Untuk 1Sis JP (syarat batas Dirichlet) 
D(i) —2x/(XP(i).sina — YP(i).cosa) 
Untuk 1S is JP 
CLi,j) - PO) # ORG, j) In OT(i,j) — PR (i, j)-In PT(i, j) — TR (4, j).AG,j) 


PO) #-2PR(i,j) PO) —PR(i,j) 4 TR, j) 








CZi,i) — : 2PO(j) Bit 
PR RN an PT) PRO TR Ai) 
Kemudian 


C(i,1) — CI(i,1)- C2(i,1) 
Untuk 2 « j « JP 


Clisi) - Cili,j) - C2(i,j) -C2(i,j 1) 
Cli, JP -1) - C2(i, JP) 
Cli, JP 42) ——2x 
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Selanjutnya ada pilihan untuk aerofoil lancip dan tumpul. 
Untuk i-JP--1 


Untuk aerofoil lancip (ya) —y (JP 41)- 0) 
cli,1) —1 
Untuk 2 Sj X JP danj - JP t2 
Cli,j) -0 
Cli, JP 1) —-—1 
Cli, JP --2) -0 
D(i)-0 

Untuk aerofoil tumpul (kondisi batas Dirichlet) 
D(i) — 2x(XP/i).sina — YP(i).cosa) 
Untuk 1S is JP 
CL(i,j) — PO(j) #-OR(i, j).In OT(i,j) - PR (4, j).In PT(i,j) — TR (4, j).AG, j) 


PO) #2PR(i,j)) PO) - PR(i,Y 4 TR(i, j)' 





C2(i, j) — 7 2P0 “In OTKi,j) 
— PR, 7 — TR(i,j) MPTi.i —PRGPTRG,)) 4. 
NN oa 
Kemudian 


C(i,1) — Ci1(i,1) - C2(i,1) 
Untuk 2S j X JP 
cli,j) — cili,j) - C2(i,j) -C2G,j —1) 
Cli, JP 1) -— C2(i, JP) 
Cli, JP --2) -—2x 
Untuk i-JPt2 (kondisi Kutta) 
Cli,1) —1 
Untuk 2S j X JP danj - JP 2 
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Clisj) - 
Cli, JP 4 1 —1 
Cli, JP 2) - 
D(i)-0 
Perhatikan bahwa untuk kasus ini 
x(j)) —y(j) untuk 1S js JP 1 
x() - wo untuk j- JP42 
IV. Metoda Panel Vorteks distribusi linear dengan kondisi batas Neumann 
Titik kontrol, T, dipilih berada di tengah-tengah panel. Kasus ekor lancip dan tumpul 


sama saja. 
Persamaan Atur 
JP 


Y Cii,i).xG) 5 DG) untuk 1 Sis JP 41 
j- 


dimana 

Untuk 1«i« JP 
D(i) -— 2n.sin(P(i) —a) 
Untuk 1 Xi « JP | 


oRfii) u Se) TR(i,j) 


CLG,j & 
Gn ia og kera) reg) 








AG) 








C2(i,j) 


TR(i,j) al Ot), ORfi,j) AG) 


PO) LPTGj)) POC) 


C3(i,j) -1-4- 








AG) 


PR(i,j) 4 Se) TRi,j) 
Po) (eta) POC) 








C4 j) 


TRi,i) 1 Ia |. si) .AGi,i) 


PO)” LPTGi)) PO) 
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Kemudian 
C(i,1) — Cl(i,1).cos AB(i,1) — C2(i,1).sin AB(i,1) 
Untuk 2 £ j « JP 


C(i,j) — CL(i,j).cos AP(i, j) - C2(i, j).sin ABG, j) 
- (C3Kf, j — 1).cos AB(i, j — 1) — C4(i,j —1).sin AP(i, j —1)) 
C(i, JP #1) - -(C3(i, JP).cos AB(i, JP) - C4(i, JP).sin AP(i, JP)) 
Untuk i-JP-1 (kondisi Kutta) 
D(i) -0 
Untuk j — 1 danj - JP - 1 
C(i,j) —1 
Untuk 2 S j « JP 
C(i,j) -0 
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BAB VIII 
PERANCANGAN BENTUK AEROFOIL 
DENGAN METODA PANEL 


Dalam bab sebelumnya kita telah membahas masalah analisis aerofoil dengan metoda 
panel. Apabila bentuk aerofoil diberikan maka metoda panel dapat diterapkan untuk 
menganalisa karakteristik aerodinamika aerofoil tersebut. Dalam masalah perancangan, kita 
ingin mencari bentuk aerofoil yang memenuhi beberapa persyaratan, misalnya mempunyai 
suatu distribusi tekanan tertentu disepanjang permukaannya. 

Sebagai contoh, data masukan yang diberikan adalah distribusi Cp sebagai fungsi 
absisa XP dari titik-titik permukaan, P, yang ordinatnya, YP, tidak diketahui dan harus dicari 
atau dihitung. Data yang diberikan adalah pasangan nilai (XP(), Cp/()) untuk 1 « JS ITA. 
Karena nilai koefisien tekanan, Cp(j), diketahui maka nilai kecepatan aliran di titik P(j) dapat 
dihitung sebagai berikut 

V(@) -“1-Cp(i) (8.1) 
Dari pembahasan di bab sebelumnya kita dapat mengambil kesimpulan bahwa untuk kasus 
aliran 2-D metoda panel yang paling efisien sebaiknya menggunakan singgularitas vortex. 
Selanjutnya diketahui pula bahwa kekuatan vortex per satuan panjang itu ternyata sama 
dengan kecepatan aliran (lihat persamaan (7.21) dan (7.22)). Ini berarti bahwa distribusi Y 
sebagai fungsi absisa XP dapat dihitung dengan mudah yaitu untuk titik P(j) yang berada 
pada permukaan aerofoil 

1) -—1-Cp(j) (8.2) 

, sedangkan untuk permukaan bawah aerofoil 

1G) - V1- Cp(j) (8.3) 
Seandainya ordinat titik permukaan aerofoil, YP(j), untuk 1 « JS JTA itu diketahui maka 
metoda panel dapat diterapkan untuk menganalisa karakteristik aerodinamiknya. Sebagai 
contoh distribusi y(j) disepanjang permukaan aerofoil dapat dihitung dengan menggunakan 
metoda panel vortex derajat 2 dengan syarat batas dirichlet. Sistem persamaan simultan yang 
harus diselesaikan adalah 


Teri: D1) -21-,.,, -D(i) untuk 1K IL JTA H1 (8.4) 


ja 
dimana untuk 1 Xi X JTA 

D(i) -2x(XT(i).sina — YT(i).cosa) (8.5) 

untuk j-1 

CK) - CL) - CX) (8.6) 

untuk 2 Si JTA-1 

Cfi, )) — CX, j) - C2i,)) #C2(i,j —1) (8.7) 

untuk j-JTA 

Cfi,j) — C2Mi,j 1) (8.8) 


Koefisien pengaruh C1(i,j) dan C2(i,j) adalah 


CX, ) - PO()) “ORG, i).In OT(i, j) - PR(i, j)-In PT, j) 








“TR, )-AGj) Bu 
Cats) POWER) POP PR AR gr.) 
4 2PO(j) (8.10) 
PR TN Pi) PN TR Ai) 
2PO(j) PO) 
Persamaan (8.4) yang terakhir adalah penerapan kondisi Kutta, yaitu 
1) --y(JTA)-0 (8.11) 


Kondisi Kutta ini tentu saja harus dipenuhi oleh distribusi y(j) yang diberikan oleh 
persamaan (8.2) dan (8.3) sebagai data masukan. 

Pembahasan secara terperinci tentang koefisien pengaruh, C(i,j), dan penurunan rumus- 
rumusnya telah diberikan dalam bab 7 dan tidak perlu diulang lagi disini. Titik kontrol atau 
titik kolokasi T(i) dengan koordinat (XTfi),YT(i)) boleh ditentukan berada di tengah panel 
PO(i) maka koordinat Ti) diberikan oleh rumus berikut 1 i X JTA —1 


151 


xTti) - 1, (xp(i) 4 xp(i-1)) 
YTti) — Y (YP(i) « Y(i 1) 


Jumlah titik kolokasi yang dibutuhkan adalah JTA-1. Ini berarti bahwa titik kolokasi yang ke 
JTA tidak berada di tengah panel, tetapi ditentukan berhimpit dengan ujung panel P(JTA), 


(8.12) 


yaitu 


XT(JTA) — akn (8.13) 


YT(JTA) — YP(JTA) 
Seandainya koordinat titik-titik T adalah data yang diketahui maka koordinat P dapat 
dihitung dengan hubungan-hubungan balik berikut. 

Untuk j-JTA 


XP('j) Ma (8.14) 


YP) -YT() 
untuk JTA —12 j21 


XP(H) - 2XT() - XTj 4 4 (8.15) 


YP(4) 5 2YT())- YT(I H-1) 
Apabila titik kolokasi P ditentukan berada di ujung panel maka hubungan antara koordinat T 
dan koordinat P menjadi sangat sederhana yaitu 
Untuk 1 Xi « JTA 


21 - ad (8.16) 
YT(i) -— YP(i) 

atau hubungan baliknya, yaitu 
sa z 2 (8.17) 
YP(i) — YTMi) 


Dalam sistem persamaan (8.4), setiap persamaan harus berbeda. Dengan perkataan lain letak 
semua titik kolokasi, T(i), itu harus berbeda satu dari yang lain. Sekarang perlu diperhatikan 
bahwa untuk aerofoil berekor lancip titik P(1) itu berhimpit dengan P(JTA) maka T(1) juga 
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berhimpit dengan T(JTA) dan ini melanggar persyaratan yang dibahas sebelumnya. Ini 
berarti untuk aerofoil berekor lancip koordinat T(1) dan T(JTA) harus dibuat berbeda dari 
koordinat P(1) dan P(JTA), misalnya dengan menentukan bahwa T(1) berada di panel PO(1) 
dan T(ITA) berada di panel PO(ITA-1) tetapi berada sedikit di depan titik ekor, yaitu 


XT(1) - XP(1)—0.001/(XP(1) — xP(2)) 2 
YT(1) — YP(1)— 0.001(YP(1) — YP(2)) 2g 
XT(JTA) - XP(JTA) —0.001(XP(JTA)-— XP(JTA —1)) ng 
YT(JTA) — YP(JTA)—0.001(YP(JTA)- YP(JTA —1)) Ba 
Hubungan balik dari persamaan persamaan di atas adalah 
XP(1) - (XT(1) —0.001(Xp(2))) /o.999 
(8.20) 
YP(1) - (ra) —0.o01( xP/(2))) Jo.999 
XP(JTA) — (xr(3Ta) — 0.001 Xp(3TA —1)))/0.999 
(8.21 


YP(JTA) — (YT(JTA) —0.001( YP(JTA — 1)))/o.999 


Sekarang marilah kita bahas masalah perancangan secara terperinci. Dalam masalah 
perancangan justru harus dicari atau dihitung. Tetapi pembahasan sebelumnya menunjukkan 
bahwa metoda panel tidak dapat diterapkan apabila koordinat tiik-titik permukaan aerofoil, 
PG), tidak diketahui. untuk mengatasi masalah ini, dalam prosedur perancangan pada 
awalnya kita harus menebak bentuk aerofoil, yaitu nilai absisa XP(j) ditentukan sedangkan 
nilai koordinat YP() terpaksa ditebak. Karakteristik aerodinamika aerofoil tebakan tersebut 


dianalisa, yaitu nilai distribusi kekuatan vortex, 1 ), dihitung dengan menyelesaikan sistem 
persamaan simultan (8.4). Nilai distribusi y(j) tersebut kemudian dibandingkan dengan nilai 
sesungguhnya yang diberikan oleh rumus (8.2) dan (8.3). Apabila nilai y(j) tebakan berbeda 


dari nilai y(j) yang ditentukan maka bentuk aerofoil tebakan, atau nilai ordinat tebakan 


YP(j) harus diperbaiki. Bentuk aerofoil tebakan yang diperbaiki ini kemudian dianalisa lagi 
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dan nilai y(j) tebakan baru yang diperoleh sebagai solusi sistem persamaa (8.4) kemudian 
dibandingkan dengan nilai y(j) yang sesungguhnya. Prosedur tersebut di atas diulang 
berkali-kali atau diiterasikan sampai akhirnya konvergen, yaitu nilai y(j) tebakan menjadi 


sama dengan nilai y(j) yang diberikan sebagai data masukan. Masalah utama dalam 
prosedur perancangan adalah mencari cara terbaik untuk memperbaiki bentuk aerofoil 
tebakan yang menjamin bahwa prosedur iterasi yang diterapkan pasti konvergen dan 
membutuhkan jumlah iterasi yang sesedikit mungkin. Untuk menentukan prosedur iterasi 
yang sebaiknya diterapkan, kita perlu melakukan analisis sebagai berikut. Pada awalnya nilai 
ordinat YP(j) yang tidak diketahui ditebak bernilai YP, (j). Nilai koefisien pengaruh CW) 
dapat dihitung sebagai fugsi koordinat (XP(j), YP,(j)) dan (XTC), YT,(j)), dimana 
YT,(j) dihitung menggunakan rumus (8.12) dan (8.13) atau rumus (8.16), (8.18) dan 
(8.19). Rumus (8.4) dan (8.5) dapat digabungkan dan diolah lebih lanjut menjadi bentuk 
berikut, untuk 1Xi X JTA. 


YT(i) -XT(i).tana - 3 Ci, 4). (9) #21. We. (8.22) 
ja 


Karena w,.,, tidak berubah dengan nilai indeks i, maka bentuk aerofoil tidak dipengaruhi 
oleh nilai w,..,- Disamping itu koefisien C(i,j) hanya ditentukan oleh bentuk aerofoil jadi 
nilai C(i,j) juga tidak dipengaruhi oleh w,.,,- Nilai w,.,, yang berbeda hanya berdampak 
menggeser letak aerofoil di arah sumbu-y. Penalaran tersebut memberikan kayakinan bahwa 
Waero boleh ditentukan bernilai sembarang tanpa mempengaruhi bentuk aerofoil yang 
dihitung. Untuk memudahkan pembahasan dan hitungan, ditentukan bahwa w,,., bernilai nol 
dan persamaan (8.22) dapat disederhanakan menjadi bentuk berikut 


YT() — XTi).tana- Y' CK,). YM) untuk 1Si  JTA (8.23) 


ja 
Seandainya nilai C(i,j) diketahui maka rumus di atas dapat digunakan untuk menghitung nilai 
YT(i) yang diinginkan secara langsung. Tetapi karena C(i,j) tidak diketahui dan dapat 
ditebak nilainya oleh C,(i,j) maka nilai YT(i) hanya dapat ditebak oleh YT, (i), dimana 
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YT, (i) - XT(i).tana - 3 C,(ij).y()) untuk 1xi« JTA (8.24) 


ja 
Karena YP,(i) ditentukan sesuka hati, sedangkan YT,(i) dihitung dengan rumus (8.24) 
yang merupakan taksiran dari solusi eksak (8.23) maka YT, (i) adalah taksiran Y'T(i) yang 
lebih baik dibandingkan dengan YT,(i). Nilai YP,(i) kemudian dihitung menggunakan 
rumus-rumus (8.14) dan (8.15) atau rumus (8.17), (8.20) dan (8.21). Nilai koefisien 
pengaruh C,(i,j) kemudian dihitung sebagai fungsi koordinat (XP(j), XP, (j)) dan 


(XT), YT,(j)). Karena YT,(i) adalah taksiran YT(i) yang lebih baik dibandingkan 
dengan YT,(i), maka C,(i,j) juga merupakan taksiran C(i,j) yang lebih baik daripada 
C,(ij). Rumus (8.23) sekarang dapat diterapkan sekali lagi untuk memperoleh taksiran 
YT(i) yang lebih baik daripada YT, (i) dengan menggunakan C,(i,j), bukannya C,(i, j) 
,sebagai taksiran C(i,j), yaitu 


YT, (i) - XT(i).tana Sad, DG) untuk 1 is JTA (8.25) 
j-l 


Prosedur di atas dapat diulang berkali-kali dengan rumus iterasi berbentuk umum sebagai 
berikut 
Untuk n— 1,2,3,.. . dan seterusnya 

JTA 


YT, (i) - XT(i).tana— Y' C, ,(ij).y(j) untuk 1 Ki JTA (8.26) 
jA 


dengan catatan bahwa nilai C,(i,j) diberikan sebagai tebakan awal seperti telah dijelaskan 
sebelumnya. 
Rumus (8.26) adalah rumus iterasi yang dapat diterapkan untuk memperoleh tebakan nilai 
YTG). Untuk nilai n yang cukup besar dapat diperkirakan bahwa 

(XT, 4 (D) — YT, (i)| Se untuk semua nilai i (8.27) 
dimana & adalah suatu bilangan positif kecil yang nilainya ditentukan sesuka hati, misalnya 
10“, dan merupakan ukuran dari ketepatan hasil iterasi atau keakuratan nilai koordinat 
YT(H) yang dihitung. Apabila n digantikan oleh (nt-1) maka rumus (8.26) dapat ditulis ulang 
menjadi 
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